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DIGITAL TELEVISION: THE USE OF WAVEFORM ESTIMATES IN ERROR 



C.K.P. Clarke, B.Sc.(Eng.), A.C.G.I. 



In p. cm. television, the video signal often consists of 
a series of 8-bit codewords which describe a "sampled-and- 
held" and quantised version of the analogue signal. If 
one of the bits in a word is altered, perhaps by distortion 
or interference, then the sample value is changed by 2^^ 
quantising levels, where N is an integer which depends on 
the significance of the bit affected. Such errors appear 
in the picture as bright or dark specks. 

Several methods have been developed to locate and 
correct errors in digital systems and these involve the 
addition of redundant information to the signal in the 
form of parity check bits. The parity bits are related to 
the modulo-2 sum (exclusive-OR function) of groups of 
bits m the signal. Hence, if an error is subsequently 
introduced, its presence can be detected by checking the 
relationships between the received parity and data bits. 
By using a system of interlinked parity checks, known as 
an error-correcting code, the error can be traced to a 
particular bit and corrected. Error-correcting codes can 
significantly reduce the impairment caused by errors in 
p.c.m. television signals,^ '^ although the level of protection 
provided depends on the proportion of parity bits to 
signal bits. 

Clearly, the addition of parity bits increases the data 
rate of the signal, so it is desirable to keep the proportion 
of parity bits to signal bits low whilst maintaining an 
adequate level of error protection. This can be achieved 
by making use of redundancy inherent in the p.c.m. signal 
and relying less on artificial redundancy added in the form 
of parity bits. 

Two kinds of redundancy are present in p.c.m. 
television signals. First there is redundancy caused by 
using a sampling frequency higher than twice the bandwidth 
allocated to the signal; for example, with a bandwidth 
of 5-5 MHz, a sampling frequency of 13-3 MHz is 
sometimes used. This over-sampling is necessary to allow 
a practicable rate of cut in the filters used in the conversions 
between the analogue and digital forms of signals; this 
redundancy is present regardless of the signal content. 
The other kind of redundancy arises because television 
signals do not always use the full bandwidth allocated 
to them. Television pictures frequently contain plain 
areas, so that most of the signal energy is concentrated at 
low frequencies and at the colour subcarrier frequency. As 
this is not true for all pictures, this second kind of re- 
dundancy is not always present. 

With sampling frequencies well above the Nyquist 
rate (for example, three-times colour subcarrier frequency) 
either or both kinds of redundancy can be exploited. 
However, with sub-Nyquist sampling at, say, twice colour 
subcarrier frequency, there is no redundancy due to over- 



sampling; in this case only the second kind of redundancy 
can be exploited. For most of this report it will be 
assumed that the sampling frequency approximates to 
three-times colour subcarrier frequency. 

Signal redundancy has already been used in error 
concealment systems-^-^ in order to reduce the effect of 
errors. !n these systems, the approximate position of an 
error is located by a simple parity check, so that the word, 
or possibly the group of words, containing the error is 
known. This word (or group of words) is then replaced 
by an estimate of the correct value, produced by assuming 
certain relationships between neighbouring sample values. 
Usually, because of the signal redundancy, these relation- 
ships hold, so that the estimate is close to the correct 
value of the word and the effect of the error is concealed. 

However, signal redundancy could be exploited by a 
system for locating and correcting errors, working either 
in combination with parity checks or, perhaps, completely 
without added parity information. The redundancy would 
be used to enable an estimate of the correct value of a 
sample to be deduced from those of neighbouring samples. 
The estimate would then be subtracted from the received 
value of the sample to produce a difference signal. If 
the estimates were reasonably accurate and no errors had 
occurred, the difference signal would be small; however, 
if errors were present, each would produce a transient 
response in the difference signal, proportional to the 
magnitude of the error. Errors could thus be located 
and corrected by extracting their positions and magnitudes 
from the difference signal. 

In essence, the production of estimates from neigh- 
bouring samples is a filtering operation, so the amplitude/ 
frequency characteristic of an estimator is a useful guide 
to its effectiveness. Therefore, in the following section, 
several methods of producing estimates are considered in 
such terms. First, the theory of estimating filters is 
outlined and this is followed by descriptions of practical 
estimators. These range from the simple filters that have 
already been used in error concealment systems to more 
elaborate filters, capable of producing near-ideal responses. 
Then, in Section 3, methods of using the estimates in 
error correction systems are described. Finally, In 

Section 4, practical applications of the correction methods 
are assessed. 



2. Sample estimation 

2.1. Theory of estimating filters 

Effective concealment or correction of an error 
depends on accurate estimation of the signal value. In 
concealment, if an inaccurate estimate were used to 
replace a word containing an error, a visible error would 
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still remain which, in some cases, could be larger than the 
original error. Similarly, if poor estimates were used in 
correction, this could lead to the wrong bit being identified 
as containing the error and would sometimes make errors 
more serious. 

Accurate estimates can be formed by passing the 
received signal through a filter with an amplitude/frequency 
characteristic which is unity at all frequencies present in 
the video signal. However, for reasons explained below, 
unit response in an estimator cannot be maintained at all 
frequencies at once, so that the amplitude/frequency 
characteristic must deviate from unity in unused parts of 
the spectrum. Thus, the presence of redundancy in the 
signal, in the form of unused parts of the spectrum, is 
essential to the formation of estimates. 



The frequency response, F(tj), of the filter is the Fourier 
transform of the impulse response, that is: 



F(u) = m e'J^t dt 



S \a,?>{t+rT) + a_^m-rT)\ e-^J^^ dt 






Estimates can be produced conveniently using a 
transversal filter in which weighted contributions from 
nearby samples are added together. The samples used 
may be taken from the line including the sample to be 
estimated, or from other lines in the same field or, 
indeed, from other fields, if sufficient storage is available. 
Therefore, one-, two- or three-dimensional estimators with 
axes of horizontal, vertical and temporal frequency may be 
contemplated. Filtering in the horizontal direction, in 
which the estimate is made from samples in the same 
line, has the advantages that the amount of storage 
required is minimal and that the horizontal bandwidth is 
usually less than half the sampling frequency. 



If the samples s_„, s„(^_ij 



>-l 



and 



5^, ^2, . . . . , s„ shown in Figure 1 are used, then Eq, 
the estimate of the value of the current sample Sq, is 
given by: 



^0 = 



n 

S 

r-1 



h^r 



-!- a 



where a_fj, . . . . , a„ are the weighting coefficients of the 
filter, as shown in Figure 2{a). Such a filter would have 
the impulse response shown in Figure 2(b), which corres- 
ponds to the equation: 



where T represents the sampling interval and b{t - rT) 
represents an impulse of unit weight atf = /"rdefined by: 



a_,j.b{t-rT) dt - a^ 



-)^-^ 



S_n ^-{n-^; 



- -^ — H — X — X— ^- 
-2 -'1 '^0 "1 2 






Fig. 1 - The relative positions of samples in a television 

line 



= S \[a,.^■a_^) cosfcoT + j(a^-a_j.) s'mrojf] 

The amplitude/frequency characteristic \F(co)\ can be 
obtained by taking the modulus of this function, and the 
phase/frequency characteristic arg[F(oj}] by taking the 

argument: 



\F(c 



n 
r=\ 



CO&-0JT 



n 
r='[ 



arg[F((u)] ^ sin 



-1 



sinrajT 



S {a^-a_j.) sinrojT 
r=1 



(i; 



(2) 



\F(co]\ 

It is convenient to make the filter symmetrical by setting: 

a_^ = fl^ ^<r<n 

as this ensures zero phase shift and halves the number of 
multiplications. So, for a symmetrical filter: 



and 



tF(oj)l = 2 S aj.cosrojT 
r=1 



arg[F(co)] = 



{3} 



(4) 



Since this expression for the amplitude/frequency 
characteristic of an estimator is made up entirely of 
cosine terms, the average level of the response curve must 
be zero; that is, the areas between the curve and the 
horizontal axis, above and below the axis, must be equal. 
Apart from this constraint, it is possible to obtain any 
shape of response, provided that the number of contribu- 
tions used is sufficient. 
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Fig. 2(a) - A transversal filter for use in estimation 



I(t) 



On O/7--I 



ai a^ 



a.^ 0_2 o_[n-^) ^-n 



-nT-{n-AM ^2T -T T 21 



[n-^)T nT f 



Fig. 2(b) - The impulse response of the filter. Note 

that the sample s^ is not used in the production of the 

estimate and that the impulse at t = rT is weighted by tiie 

coefficient a ^ 



The ideal amplitude/frequencv characteristic for an 
estimator would be a response of unity from to 5-5 MHz, 
as this would ensure that all samples in legitimate video 
signals could be accurately estimated. To satisfy the 
constraint on the average level, the characteristic must 
reach large negative values in the out-of-band region 
between 5'5 IVlHz and half sampling frequency (6-65 IVlHz 
for three-times subcarrier sampling), as shown in Figure 3. 
In theory, the shape of the characteristic in the out-of-band 
region is not important although, in practice, a curve 
which falls away gradually to a large negative value at the 
half sampling frequency would be most suitable. This is 
iDecause the signal may still contain components slightly 
above 5-5 IVlHz, as the pre-sampling filter in the analogue- 
to-digital converter still has a finite response in the out-of- 
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band region. Further, a gradually falling characteristic may 
be obtained with fewer stages than a more abrupt cut-off. 

The characteristic of Figure 3 only makes use of the 
signal redundancy due to over-sampling. However, as 
mentioned in the introduction, video signals often contain 
more redundancy than this because, for much of the time, 
the signals are confined to low frequencies and to the 
region around colour subcarrier frequency. For signals 
that conform to this frequency distribution, a filter with 
the response of Figure 4 will produce accurate estimates. 

2.2. Estimators used in error concealment systems 

Characteristics approximating that of Figure 4 can 
be produced using very simple filters and some of these 
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Fig. 3 - The ideal amplitude/frequency characteristic for 
an estimator 



Fig. 4 - An estimator characteristic suitable for signals 
with no significant high-frequency luminance content 
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have already been used in error concealment systems. For 
example, in one 8-bit p.c.m. system the value of the 
third-previous sample was used as the estimate.-^ This 
corresponds to an asymmetrical transversal fitter with a 
single non-zero coefficient, ii_2, equal to one. For this 
filter, from equations (1) and (2); 



\F 



CO 



= 1 



and 



arg[F(u))]= -3coT 



Although the amplitude response is unity, the estimates 
are accurate only when the phase response is zero or an 
integer multiple of In radians. 
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The effect of the phase shift is more apparent when 
the value, as a function of frequency, of the difference 
between the estimate and the current sample is considered; 
this characteristic, shown in Figure 5, is calculated from 
the following expression: 

\F{co]\ = (2-2cos3wr)'''= 

vi^iich can be obtained by a derivation similar to that for 
equation (1). With a sampling rate of three-times 

subcarrier frequency the estimator gives perfect accuracy 
(zero error) at zero frequency and at subcarrier frequency. 
However, the amount that the characteristic deviates 
from zero shows, for sinusoidal signals, the magnitude 
of the maximum error that can occur between the 
estimate and the actual value. 

The generally poor accuracy of this estimator limited 
the effectiveness of the concealment. Therefore, the 
parity bit used to protect each word was formed using 
only the three most-significant bits of the word. Errors 
in the five least-significant bits, since they were not 
included in the parity check, were never detected and 
so were not concealed. Applying concealment to errors 
in less-significant bits was found to increase the impairment, 
because the estimates were generally less accurate than 
words containing low-significance errors. Even so, this 
system provided a useful reduction in the level of 
impairment caused by errors. 
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Fig. 5 - The magnitude of the maximum error produced 
by a ttiird-previous-sample estimator 



Fig. 6 - Tiie amplitude/frequency characteristic of the 

average of two samples, each spaced one subcarrier period 

away from the sample under consideration 



Another error concealment system used an improved 
estimator which averaged the samples one subcarrier 
period before and after the current sample. Although 
two samples were used in the calculation, this produced 
a symmetrical filter in which only one coefficient value, 
and therefore one multiplication, was required. 



The amplitude/frequency characteristic of this estima- 
tor, which was used with sub-Nyquist sampling at twice 
subcarrier frequency, is shown in Figure 6. This 

characteristic is obtained by substituting ^2 ^ ^/^ '"^ 
equation (3) to produce: 



\F{(jj)\ = co%2coT 

In order to obtain the same characteristic with 
sampling at three-times subcarrier frequency ^3 = V2 should 
be used instead. For symmetrical filters there is no 
relative phase shift between the estimate and the sample 
under consideration and, therefore, the accuracy of the 
estimator is determined by the closeness to unity of the 
amplitude/frequency characteristic. Therefore, this esti- 
mator provides a genera! improvement in accuracy over 
the third-previous-sample estimator, especially at low 
frequencies and near subcarrier frequency. Consequently, 
this estimator allowed the parity protection to be extended 
to cover the first four bits of a word, counting from the 
most significant bit. 

2.3. Improved estimators 

A worthwhile improvement in response over the 
estimators described in the previous section can be obtained 
from reasonably simple symmetrical filters. For example. 
Figure 7 shows a characteristic, calculated from equation 
(3), that can be produced using four coefficients. This 
estimator is much more accurate in the range to 
3-75 MHz than that of Figure 6 and, although it is less 
accurate for the subcarrier sidebands, unity response at 
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Fig. 7 - The amplitude/frequency characteristic of an 

estimator using four samples on either side of the sample 

under consideration. Coefficient values: 

a^ = 0-940, fl2 = -0-766, flg = 0-500 and ^4 - -0-174 

subcarrier frequency has been retained. 

An alternative characteristic, again obtained using 
four coefficients, is shown in Figure 8. Here, some 
accuracy at low frequencies has been sacrificed in order 
to use simple coefficients, with which multiplication can 
be accomplished merely by shifting the significance of 
the signal bits. 

Further improvements in the estimator characteristic 
can be gained by increasing the number of samples 
gontributing to the estimate. If sufficient contributions 
are used, characteristics similar to the ideal of Figure 3 
can be produced. Until recently, the complexity of 
these near-ideal filters would have prevented their use; 
however, because read-only memories can now replace 
the arrays of adders previously used for multiplication, 
tfie production of such filters is relatively simple. 



describing a general polynomial of order 2h-1 passes 
through 2n known sample values (n on either side of the 
sample under consideration). Then the coefficients a^., 
r samples are away on either side of the sample under 
consideration are given by: 



(-1 



v^\ 



[n\ 



a^ ^ 



(n+r]\ {n-r)] 



The derivation of this formula is given in the Appendix, 
Section 7.1. The amplitude/frequency characteristic 

given by these coefficients, shown in Figure 9, has the 
merit of being flat in the pass-band; however, a high 
value of n is required to maintain a reasonable response 
to 5-5 MHz. For high values of r, the coefficients are 
very small and can be assumed to be zero without 
noticeably affecting the shape of the characteristic. For 
example, for n = 40, setting a^-j to a^Q to zero has a 
negligible effect on the characteristic, so that only sixteen 
coefficients are needed. 
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Fig. 9 - A set of amplitude/frequency characteristics pro- 
duced by curve-fitting using a polynomial of order 2n— 1 

with different values of n 



Three methods of calculating the weighting co- 
efficients have been investigated, each yielding a different 
approximation to the ideal characteristic. The first method 
consists of calculating the coefficients so that the curve 
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Fig. 8 - The amplitude/frequency characteristic of an 

estimator with four coefficients chosen for ease of 

multiplication: 

fli = 1,fl2 ^ ^^'^''^3 ^ y3anda4 = —% 



(a) /( - 10 
Ic} n = 40 



(b) n = 20 
(d) n - 80 



For the second method of calculating the filter 
coefficients, the ideal characteristic was assumed to be 
flat in the pass-band and then to follow a cosine law to a 
large negative value at the half-sampling frequency. The 
coefficients for this shape of characteristic are calculated 
from: 



sin 2-nrf T 



a, = 



where /„ represents the width of the pass-band and T again 
represents the sampling interval. The derivation of this 
formula is given in the Appendix, Section 7.2. In general, 
an infinite number of coefficients is needed to achieve the 
ideal characteristic, but those for the higher values of r 
can be neglected. With sixteen coefficients, this method 
gives an extended high-frequency response, although the 
characteristic contains ripples in the pass-band. This is 
shown in Figure 10. 
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The main disadvantage of such a concealment system 
is the allocation of one parity bit per word, which is a 
higher proportion than is needed for some error correcting 
codes. Furthermore, estimates can be used more 

effectively as a means of locating errors, thus reducing 
the requirement for parity. This is described in the 
following section. 



3.1. Correction methods using simple estimators 



Fig. 10 - The amplitude/frequency characteristic obtained 
by using the first 16 coefficients of a filter calculated to 
produce a cosine law in the region above 5-5 MHz 

For optimum performance using television signals, the 
accuracy, of the estimator should be related to the power 
spectrum of the signals. In general, this spectrum is not 
well defined, but a worthwhile advantage would be gained 
by ensuring better accuracy at low frequencies and 
around the subcarrier frequency than at less important 
signal frequencies. Therefore, a third method of calculating 
the coefficients was derived, which allows the response 
values to be specified at particular frequencies. As the 
response value at any frequency w is given by: 






{3} 



it is possible to write n such equations and solve them for 
the n unknowns a^ to a^. Thus the response value \F{co)\ 
can be set at n frequencies co^ to oj^j when an n- 
stage fitter is used. In practice, the frequencies used must 
be approximately evenly spaced as, otherwise, the response 
actually obtained may deviate wildly from the required 
response in regions where the set points are widely 
separated. However, it is possible to obtain solutions 
that represent a compromise between the responses of the 
other two methods and to ensure that the response is 
accurate at the more important frequencies. Figure 11 
shows a characteristic obtained by this method using 
sixteen coefficients. 



The estimating fitters described in Section 2 can be 
used with advantage for the detection and location of 
errors. Thus, by mat<ing efficient use of the redundancy 
inherent in the signal, the proportion of parity bits 
required can be reduced. 

Although the best performance would be obtained 
with near-ideat estimators, the simple estimators already 
used for concealment could be used to locate and, 
subsequently, correct errors. This would eliminate most 
of the residual impairments caused by using inaccurate 
values for concealment, because error location is less 
dependent than concealment on the accuracy of the 
estimates. Therefore, a worthwhile improvement in 
performance could be gained. However, comparison of 
the estimates with the received values does not, by itself, 
give a reliable indication of the presence of errors. For 
this reason, a parity checking system must be used in 
addition, to ensure that only words likely to contain 
errors are altered. 

Two methods of error correction based on simple 
estimators have been devised and these use different 
parity systems. The methods will be explained, and 
compared with concealment, in a series of worked 
examples using the signal of Figure 12. In each case, 
the estimates will be formed by averaging two samples 
occurring three sampling intervals before and three after 
the sample under consideration (the amplitude/frequency 
characteristic of this process is shown in Figure 6). This 
estimator is unlikely to give accurate values because the 
signal in Figure 1 2 has a large component at approximatelv 



2.4. Error conceatment using accurate estimates 

If the near-ideat estimators described in the previous 
section were used in an error concealment system, almost 
ail the errors could be successfully concealed. This 
would allow the parity protection to be extended, perhaps 
to include all the bits of each word. However, with the 
simple parity system normally used in error concealment, 
two mechanisms would remain through which errors could 
affect the signal. First, if two errors were to occur in 
neighbouring words, each would cause an error in the 
estimate used to conceal the other. Secondly, assuming 
simple parity, pairs of errors would not be detected if they 
occurred in the same word, or if there was one error in a 
word and another error in the associated parity bit. Even 
so, with random errors, these imperfections would only 
become noticeable at high error rates. 
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Fig. 1 1 ' The amplitude/frequency characteristic for a set 

of coefficients calculated to fix the response value at 

unity at sixteen equi-spaced frequencies ranging from 

to 5-5 MHz 



(PH-178) 



-6- 




J L_L 



J L_i. 



M-^ M-4 M-Z M M^Z A^t4 /l^+6 

{^ 

sample numbers time 



for the incorrect value (11) and this would leave an error 
of 72 levels. This error, caused by the poor accuracy of 
the estimate, is larger than the original error of 64 levels. 

Alternatively, in the correction method, the estimate 
is used in an attempt to select the original value of the 
vi'ord, instead of substituting the estimate directly into the 
signal. The method makes the assumption that, when a 
parity check fails, only one error has occurred in the 
protected bits. This is a reasonable assumption for random 
errors, because even numbers of errors will not be detected 
and odd numbers of errors greater than one virill occur 
extremely infrequently. Therefore, the correct word will 
be among the five words obtained by altering the protected 
digits of the received parity group one at a time. In the 
example, the five possible values for word Af are: 



Fig. 12 - A section of sampled video waveform tal<en from 
the diagonal bars of a test card F signal 

1-5 MHz. it will be assumed, for the examples, that an 
error has occurred in bit 2 of word M, changing its value 
from 75 (01001011) to 11 (00001011). 



3.1.1. A correction 
each word 



method using one parity bit for 



With a parity system in which one parity bit is 
formed from the four most-significant bits of each word, 
the estimates can be used either for correction or for 
concealment. In the example, the error in word M is 
detected by the following parity check on the first four 
bits of word M: 



or 



139 (10001011) if the error is in bit 1 

75 (01001011) if theerrorisin bit2 

43 (001 01 01 1 ) if the error is in bit 3 

27 (00011011) if theerrorisin bit4 

11 (00001011) if the error is in the parity bit 



The value that is most likely to be correct is chosen from 
these five by selecting the value nearest to the estimate. 
In this case, the value of the estimate is 147 so that 139 
would be chosen. Although this is slightly closer to the 
correct value than was the value substituted by the 
concealment method, an error of 64 levels still remains. 

3.1.2. Re-assigning the parity bits 



o®oeo®o ^ 1 

where ffl signifies modulo-2 addition. The bits involved in 
this parity check are shown ringed in Table 1. Next, an 
estimate 

TABLE 1 

sample decimal binary value parity 

number value 2^ 2^ 2^ 2^ 2^ 2^ 2^ 2° bit 



10 10 1 

110 10 1 

10 10 10 1 

® ® © © 1 1 1 ® 

110 10 11 

10 10 1 1 

10 11110 



of the value of word M is formed by averaging words 
M~3andM+3: 



M-2 


136 


M-2 


105 


M~1 


74 


M 


11 


M+-\ 


107 


M+2 


137 


M+3 


158 



An improved correction method can be produced 
by abandoning the allocation of one parity bit to each 
word. Instead, an arrangement can be used in which each 
parity bit is formed from a group of bits of the same 
significance from several words.^ The parity groups are 
interleaved so that the two words used to produce estimates 
{occurring three sampling intervals before and three after 
the sample under consideration) are always taken from 
groups different from the one in which an error has been 
detected. This ensures that the error does not confuse 
the selection process by producing wrong estimates for 
words belonging to the same parity group. As before, 
only the four most-significant bits of each word are 
protected and this produces the same ratio of parity bits 
to data bits as for the previous parity system. The 
parity bits are calculated according to the following rules, 
with the notation defined in Table 2; 






y. (136-1- 158) = 147 
In a concealment system, the estimate would be substituted 
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P^y ^ b^ + d^ + f^ + h, 



TABLE 2 

sample decimal binary value parity 

number value l^ 2® 2^ 2^ 2^ 2^ 2^ 2^ bits 



averaging the words three sampling intervals avi'ay on 
either side; 



2 


h 


3 


c 


4 


d 


5 


e 


6 


f 


7 


g 


8 


h 



«1 «2 '^S ^A ^5 «6 '^y «£ 
^1 ^2 *3 ^4 ^5 ^6 h h 



t?1 (^2 ^3 c/4 ^5 t^6 "^7 "^8 

^1 63 £3 64 £5 eg ey eg 

/1 /2 /s A -^5 -^6 fl fs 

^1 ^2 ^3 S4 S5 Sq Si Sz 

h-^ hj ^3 h^ //5 Hq h-i h^ 



P^x 

P2x 
P2y 
P3X 
P3y 
Pax 
PAy 



Applying this parity system to the signal values used 
in the example produces the parity bits shown in Table 3. 
The presence of an error is indicated by the p2y parity 
check, since 

1 fflO®0®0 i 

The appropriate bits are shown ringed in Table 3. 

TABLE 3 



sample decimal 
number value 



Af-6 
M-5 

M-4 



161 
165 
161 



binary value parity 

27 26 25 2^ 23 22 2l 20 bits 

10 10 1 

10 10 10 1 

10 10 1 



M-3 136 10001000 

M-2 105 (D 1 1 1 

M-\ 74 10 10 10 

M 1,1 ® 1 1 1 © 

M+1 107 110 10 11 1 

M+2 137 1 (0) 1 1 

M+3 158 10 11110 1 

M+4 167 1 (0) 1 1 1 1 

M+5 165 10 10 10 1 

M+6 164 10 10 10 

M+7 163 10 10 11 

The parity check shows that the error could be in the 
second bit of any of the four words M-2, M, M+2, M+4. 
Therefore, two possible values should be considered for 
each of these words and these are; 

105or 41 (01101001 orOOIOIOOl) forwordM-2 

11 or 75 (00001011 or 01001011) forwordM 

137 or 201 (10001001 or 11001001) for word M+2 

167or231 (101001 1 1 or 1 1 1001 1 1) forwordM+4 

Then, an estimate is formed for each of these words by 
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Va (165+107) 

V2 (136 + 158) 

Va ( 74 + 165) 

V2 (107+163) 



136 for word M-2 
147 for word M 
11 972 for word M+2 
135 for word M+4 



The two possible values for each word are compared with 
the estimate and the nearer of the two selected for 
substitution into the signal. In this case, the values 
105, 75, 137 and 167 are selected and this eliminates 
the error in word M. If in alt four cases the estimates 
had been closer to the received values of the words than to 
tiie alternatives, it would have been reasonable to suppose 
tfiat the error had occurred in the parity bit and therefore 
to take no action. 

Although, in the example, this method has corrected 
the error, the method will not always be able to work 
with very inaccurate estimates. However, in general, the 
method is likely to be a considerable improvement over 
the concealment and correction methods described above 
using one parity bit with each word. The strength of the 
method is that the probability of correcting an error 
increases with the significance of the bit containing the 
error. Thus, for errors in the most-significant bit, which 
are very noticeable, the probability of correction is very 
high. 

In general terms, the estimates will have a probability 
mass function centred on the actual value of the word; 
this is illustrated in Figure 13. The mass function will 
vary both with the value of the word and with changes of 
picture content. If the word is altered by an error, the 
received value will be displaced from the actual value by 2 
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Fig. 13 - The relative positions of ihe actual value, the 
received value and the probability mass function of the 
estimates on a horizontal scale of sample values. The 
shaded region shows the proportion of estimates that are 
not sufficiently accurate to correct the error 



levels. The error will be corrected only if the estimate is 
nearer to the actual value than to the received value. 
Therefore the probability that the error will be corrected 
Is the sum of all the probabilities under the probability 
mass function in the unshaded region. Clearly, the 
further the received value is displaced from the actual 
value, the greater is the probability that the error will be 
corrected because, in such cases, estimates that are more 
inaccurate can be tolerated. 

Another advantage of the method is that the 
proportion of parity bits used is not tied to the equivalent 
of one bit per word. A smaller ratio of parity bits to 
data bits could be used by increasing the number of data 
bits in each parity group. Naturally, this would tend 
to degrade performance but, for situations in which the 
error rate is suitably low, this system could afford 
adequate protection with a very small proportion of 
parity bits. 

However, these methods of correction using simple 
estimators were considered before the possibility of 
producing near-ideal estimators had been realised. It 
now appears that, if estimates of consistently high accuracy 
can be formed, the selection processes involved in correction 
will become relatively unimportant. Even so, simple 
estimators could be used with selection when sub-Nyquist 
sampling is used. 

3.2. Correction methods using near-ideal estimators 

3.2.1. Error location and correction 

The highly accurate estimators described in Section 
2.3. could provide a means of correcting errors without 
the need for parity bits. A reliable indication of the 
presence of errors can be obtained by subtracting the 
estimates from the received signal values to produce a 
difference signal. Then, if the estimates and received 
values are not in close agreement, a large difference is 
produced which indicates that the received signal contains 
an error. 

The combination of the estimator and the subtractor 
forms a high-pass filter which, ideally, gives no output for 
video signals; Figure 14 shows the high-pass characteristic 
obtained by combining the estimator of Figure 3 with 
a subtractor. The occurrence of an error in the p. cm. 
signal is equivalent to the addition to the video waveform 
of an impulse function, which contains energy at all 
frequencies. Therefore, the filter, by allowing through the 
high-frequency content of the impulse and rejecting the 
components of legitimate video signals, produces an output 
only when an error is present. 

The difference signal resulting from an error in the 
video signal has the same form as the high-pass filter 
impulse response, scaled according to the significance of 
the bit in error. The precise form of this signal depends 
on the type of estimator used, but all near-ideal estimators 
will produce waveforms of similar shape. For example, 
v^en the estimator of Figure 11 is used, the waveform 
produced by an error will have the same shape as the 
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Fig. 14 - An ideal fiigh-pass cttaracteristic produced by 
subtracting the estimate from the received value 

impulse response shown in Figure 15. 

The difference signal waveform contains all the 
information needed to locate the bit in error. The 
central impulse marks the position of the word containing 

the error and the magnitude of the central impulse is 
equal to the significance of the bit in error. The sense 
of the error is also given: when a has been changed to 
a 1, the central impulse is positive, whereas for a 1 to 
error, the central impulse is negative. Therefore, the 
error location process consists of extracting this information 
from the difference signal waveform. 

A simple method of locating errors would consist of 
detecting peaks in the magnitude of the difference signal 
(that is, the difference signal with the sign suppressed). 
For example. Figure 16 shows the "magnitude waveform" 
produced by three errors, two of which have occurred 
in the same word. The height of the first peak, 32 levels, 
shows that an error has occurred in the third-most- 
significant bit of the corresponding 8-bit word, whilst 
the second peak, which has a height of 80 levels, shows that 
there are two errors, one in the second bit which has 
contributed 64 levels and the other in the fourth bit 
which has added 16 levels. When the position of an error 
has been found (the peak of the magnitude waveform) the 
sense of the error can be obtained directly from the sign 
of the central impulse. Finally, to correct the error, an 
impulse of equal magnitude and opposite sign would be 
added to the received value of the word. 

3.2.2. Improving noise immunity in the Eocation 
process 

In practice, the difference signal might be degraded 
by noise. This would be produced by video and noise 
components above 5-5 MHz in the video signal, by 
quantising errors and by legitimate video information 
that has not been fully rejected. Although components 
above 5-5 MHz are not expected to be a problem, if 
significant out-of-band components are present their am- 
plitudes could be reduced by a digital low-pass filter 
immediately following the analogue-to-digital converter. 
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Fig. 15 - The impulse response formed when the estimator of Fig. 1 1 is combined with a subtractor 



Even so, some noise will be present and this could affect 
the performance of the peak detector method. 

A better method of finding the central impulse 
would be to form the average of groups of m samples 
(1 < m < 16) and then to subtract the averages formed 
m + 1 samples apart. This process would have the 
impulse response shown in Figure 17 which, because it is 
anti-symmetric, produces a 90° phase shift at all frequen- 
cies. This would act as a differentiator and would produce 
an output related to the slope of the input waveform. For 
the difference-magnitude waveform produced by an error, 
the differentiator output would be zero at the central 
impulse and would have large values on either side. 
Therefore, detection of the minimum slope in the difference- 



magnitude waveform would provide a sensitive method of 
obtaining the position of a word containing an error. 
Because the process involves averaging, detection would be 
less affected by noise in the difference waveform, especially 
for large values of m. 



Similarly, the effects of noise could be reduced by 
the use of an averaging process to obtain the magnitudes 
of errors. Thus, the values at several points in the 
difference waveform could be used to find the magnitude 
of an error, rather than using the value of the central 
impulse directly. This is possible because each impulse 
in the difference-signal waveform is related to the mag- 
nitude of the error by the appropriate fitter coefficient. 
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Fig. 16 - The magnitude of the dif- 
ference signal produced by three errors, 
two of which have occurred In the 
same word 
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Fig. 17 - The impulse response of a band-pass differentia- 
tor which could be used to improve the accuracy of error 
location 

3.2.3. The effect of errors in nearby words 

As described so far, the method is not capable of 
correcting errors that have occurred in words separated 
by less than the "width" of the estimating filter. This 
is because the difference-signal waveforms produced by 
such errors would overlap and this could confuse the 
error-location logic. However, there are several ways 
in which the method could be extended in order to deal 
with these errors. 

One method would consist of providing a second 
set of estimates based on samples from other lines. Then, 
if the estimates from the horizontal filter did not clearly 
indicate the error positions, the vertical error indications, 
which would not overlap, could be used instead. The 
general use of vertical estimates has two main dis- 
advantages: first, it is costly, because several line delays 
are required to make a sharp-cut filter and, secondly, 
there is no inherent redundancy in the vertical direction, 
because vertical frequencies can exceed half the vertical 
sampling frequency corresponding to the lines on one 
field. However, for back-up estimates a relatively 

simple vertical filter would be sufficient and the lack of 
redundancy would not be a problem. 

Alternatively, a small proportion of parity bits 
could be added to the signal. Each parity bit would 
protect all the bits of several words, with the words 
separated by at least the width of the filter. Then the 
positions of the words containing the errors could be 
found and the errors partly corrected by replacing these 
words by estimates. This method has the disadvantage 
that each estimate would still be affected by the other 
error; so, for adjacent or nearly adjacent errors, the 
accuracy would be poor. 

Further investigation may well produce attractive 
alternatives to the methods described above. 



Although many of the operations involved in using 
estimates for error correction are untried, their performance 
is entirely predictable because they can be implemented 
with digital circuits. However, assumptions have been 
made about the properties of digital television signals. 



Therefore it is unreasonable to attempt quantitative 
predictions of the reductions in impairments due to errors 
that could be obtained. Even so, it seems likely that a 
worthwhile reduction in uncorrected error rates could 
be obtained entirely without parity protection. 

If parity protection is found to be unnecessary, 
no special equipment will be required at the sending 
end of a channel. Therefore, the correction system using 
estimates could be applied in situations where errors had 
not been expected to occur, or applied after error- 
correcting parity code systems to obtain a further reduction 
in error rate. 

The system could be used for the correction of 
burst errors and random errors. Short burst errors 
(less than one television line) could be dealt with by 
using horizontal estimates for error detection and vertical 
estimates for error location. Further, burst errors in one 
channel of a "de-multiplexed" system (for example, a 
multi-track tape recorder) could be corrected. This would 
be particularly straightforward for systems in which the 
number of sub-channels exceeds the number of samples 
in the width of the estimating filter; the waveforms 
produced in the difference signal by individual errors 
in a burst would be separate. For systems with fewer 
sub-channels, it would be possible to place zeros in the 
impulse response of the estimator to avoid positions 
containing other errors belonging to the same burst. In 
general, when the occurrence of errors is not completely 
random, some benefit could be gained by designing the 
correction system to take advantage of the characteristics 
of the errors. 

When used with signals sampled below the Nyquist 
rate, the error protection methods must assume that some 
parts of the video spectrum are little used. Error con- 
cealmenf^ has already been used with sub-Nyquist signals^ 
with reasonable success. However, because there is no 
out-of-band region with these signals, there is less scope 
for producing improved estimators. 

An estimator could be produced with a narrow, but 
very deep, notch at some little-used frequency, probably 
in the region of 3 MHz; however, if there was energy 
present at this frequency, its effect would be greatly 
amplified, because of the depth of the notch, so that 
the estimates produced would be of poor accuracy. 
Therefore, in sub-Nyquist systems, it would be necessary 
to make full use of the selection methods described in 
Section 3,1., and of vertical estimates, in order to reduce 
the reliance on possibly inaccurate horizontal estimates. 

Although this report has been dealing with television 
signals, the principles involved also apply to error protec- 
tion for p.c.m. sound signals. Indeed, several simple 
estimators have already been used in an error concealment 
system for sound signals.^ 

In general, an estimator for sound signals would 
require more stages than an estimator for video signals. 
This is because sound signals normally occupy a greater 
proportion of the spectrum of the digital signals (15 kHz 
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bandwidth with a half sampling frequency of 16 KHz for 
sound, compared with 5-5 MHz bandwidth with a half 
sampling frequency of 6-65 MHz for video). Therefore, 
the out-of-band region for sound is smaller, so that a 
faster rate of cut is required to achieve the ideal response. 

However, although the filtering operations involved 
in estimators for sound signals would be more complex, 
iheir realisations could be simpler than those for video 
signals; the distributed processing that would be used in 
video digital filters could be replaced by a single central 
processor performing several hundred operations in each 
sample period of a sound signal. When suitably fast 
devices become available, these operations could be 
performed by a microprocessor. 



5. Conclusions 

Methods of protection against errors in digital systems 
are usually based on the addition of redundancy to the 
signal in the form of parity check digits. An alternative 
approach has been presented which makes use of redundancy 
inherent in the digital signal, and thus avoids the increase 
in data rate required to accommodate parity digits. 

Methods based on inherent signal redundancy make 
use of estimates of each sample value produced from 
neighbouring samples, using a filter. Several ways of 
synthesising very accurate estimating filters have been 
described. It seems that filters using sixteen samples 
preceding and sixteen following the sample under con- 
sideration will produce estimates of very good accuracy. 

One application for these filters would be in im- 
proved error concealment systems, replacing the simple 
estimators that have been investigated previously. -^-^ 
However, by using the estimates for error location, 
errors could be corrected without the need for parity bits. 



Methods of selecting the most-probable values of samples 
have also been described. These could improve the 
effectiveness of the error location methods, especially 
if the accuracy of the estimates is generally poor. 
However, with estimates of consistently high accuracy, 
the selection methods may be superfluous. 
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7. 

7.1. The derivation of (»efficients for transversal filters based on polynomial curve fitting 

An estimate of tlie current sample value (Vq) can be found by summing weighted proportions of the n sample values 
on either side of the current sample: 

yo = (i^ny^n +■■■■+ «-i7_i + «1^1 + ■ ■ ■ -^"nyn 
where a_„ . . . .a__^ anda^ . . . . fl„ are the weighting coefficients. If;' is chosen to lie on a polynomial of order 2m-- 1: 

y = (x-x^) {X-X2) U-X2,j_i) 

where x^^2' - ■ ■ ■ , ^2n-^ ^""^ ^^^ 2n—^ solutions of the polynomial, this expression for;' can be substituted into the 
equation for^'Q to produce: 

(x-x-^) {X-X2) .... {x~X2n-^) =- «_n (x-n-x^) {x-n-X2) .... U-«-X2„_i) 

+ «-(o_-l) {x-n+^-x■^) {x-n+'\-X2) .... {x-n+^-X2n-^) 

+ ....+ 

+ fl„ (x+n-x.^) U+M-X2) .... U+«-X2r--i' 

The values of the coefficients fl_„ to a„ required to make the polynomial pass through the 2« known sample values can 
linen be obtained by setting x^ to X2„--\ to values which make all but one of the terms on the right-hand side zero. So 
for the rth. coefficient: 

Xt = x-n 



'■fi+1 



= JC-1 

= x+^ 



x„^. = x+r+^ 



^2n-^ 



x+n 



so that: 



{x-x+n) ... {x-x-i-^) {x-x-^) ... U-x-r+l) U-x-r-1) ... (x-x-n) 

«r = ~ — ^ 

(x+r-x+n) . . . (x+r-x+^] {x+r-x-1) . . . {x■^r-x-r+^) (x^-x-r-^) , . . (x+r-x-n) 

n ... 1 {-^)...{-r+^)(-r^■\) ... (-n) 
(«+»■)... (r+1) {r-D... 1 (-1} ...i-n+r) 

_ (-ir+i («!)2 

{n+r) I (n-r) ! 
7.2. The derivation of coefficients for transversal filters with cosinusoidal stop-band responses 
The shape of the required response is defined by: 



\F{f)\ = 1 



_ cos— ^^ 2fT 

i-2f„r L i-2f„r 



-2/pr L 1-2/p. 

(PH-178) "13- 



2T 



where/ represents the width of the pass-band and T represents the sampling interval. 
If the response shape is resolved into cosine components so that: 

\F{f)\ = 2fl^cos 2nTf + 2a2COsAirTf + .... + 2a^cos2TrrTf + 



then 



a^ = T 



»1 
2T 

-1 
2T 



\Fif)\ cos2TrrTf df 



r = ^,2, 



= 2T 



2T 



\F{f) I cos2nrTf df 



«/ 



If 2T 

= 2T I ^cos2nrTf df + 
J 



'fp 2T 



= 2T I cos,2-nrTf df + 
J 




cos 



2T^{f-fp)T 
^'-^fpT 



2fpT 



cos2-nrTf df 



2n{f-fp)T 
cos — — - — -.cos2TTrTf df 



^-2fpTjfp ^-2fpT 



cos27TrTf df 



^-2fpTjfp 



After integration this becomes: 



1 r{]-2fpT) 
— + _ — + 



2/^7- 



nr ^r{^-r^ (1-24712] M^-7fTl 



s\n2TirTf 



s\n2TrrTf„ 



TrrC\-2fT) ^-r^^^-2f T)^ 
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